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Abstract: After updating the determination of the combination of Kobayashi-Maskawa 
elements Im according to our new estimate of the parameter ]3k, we study the 

OF-violating ratio e'/e by means of hadronic matrix elements computed to O(p^) in 
the chiral expansion. It is the hrst time that this order in chiral perturbation theory is 
included. We also discuss the relevance of some O(p^) terms that are generally neglected 
in the calculation of the electroweak penguin matrix elements. The most important effect 
of this improved analysis is the substantial reduction (20%) of the leading electroweak 
penguin matrix element and accordingly a reduced cancelation between the electroweak 
and gluon penguin contributions. The ratio e'/e is thus larger than previously estimated 
and its predicted value enjoys a smaller uncertainty. We hnd e'/e = 1.7 X 10-T 
Values positive and of the order of 10“^ are therefore preferred. 
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1 Introduction 

The real part of the ratio e'/e is a fundamental parameter of the standard model that 
measures the amount of direct CP violation in the decay of a neutral kaon in two pions. 

Most recent estimates (see refs. [1, 2, 3]) have so far agreed that values of the order of 
10“^ are to be preferred within the standard model. The reason for the smallness of these 
values resides in the cancelation between the contributions of two classes of diagrams: the 
gluon and the electroweak penguins. Moreover, because each class gives a contribution 
of order 10“^ and of the opposite sign, the resulting value of order 10“^ is plagued by an 
uncertainty of at least 100%—as it can be readily understood by assigning an optimistic 
uncertainty of 10% at each of the two contributions. 

Our previous estimate of e'/e [3], was based on the leading order (LO) lagrangian 
containing 0(p°) and 0{p^) terms, where the corresponding coefficients were calculated by 
means of the chiral quark model (yQM) [4, 5]. The O(p^) mesonic one-loop corrections 
for the 7171 matrix elements were then calculated via the LO chiral lagrangian. 

Finally, the £t of the AJ = 1/2 selection rule in tttt decays allowed us to constrain 

the ranges of the free parameters involved. Nonetheless, in ref. [3] we were forced to 
assign a rather large uncertainty to our prediction of e'/e because of the almost complete 
cancelation between gluon and electroweak contributions. 
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We now take our analysis one step further by including, in the calculation of the 
hadronic matrix elements, the terms proportional to the quark current masses and up to 
O(p^) in momenta. These are terms next-to-leading order (NLO) in the AS = 1 chiral 
lagrangian determining e'/e . 

As we shall discuss, the inclusion of the terms proportional to the current quark masses 
and those of O(p^) shows that the cancelation we hud in the LO is substantially reduced. 
As a result, the preferred central value for e'/e is now above 10“^, even though values of 
e'/e of order 10“^ cannot be excluded. 

In summary, the new elements that we have introduced in the present analysis are: 

• a term proportional to the quark current masses to be added to the 0{p^) 
corrections to the leading 0(p°) part of the lagrangian which characterizes the 
electroweak penguin contribution; 

• the O(p^) matrix elements of the relevant AS = 1 four-quark operators directly 
computed in the xQM via quark integration (these matrix elements correspond to 
linear combinations of chiral coefficients of the O(p^) chiral lagrangian). This order 
has never been included before; 

• the contribution of the chromomagnetic operator Qn (which arise at 0(p^)), first 
discussed in ref. [6] but so far never systematically included; 

• updated values for asijnz) and mt] 

• updated ranges of the xQM parameters {qq) and {agGG/Ti) obtained from our 
recent O(p^) £t of the AI = 1/2 selection rules in kaon decays [7]. 

The most important consequences of these improvements are summarized by the fol¬ 
lowing results: 

• the effect of consistently including all O(p^) corrections to the 0(p°) term in the LO 
chiral lagrangian reduces by about 20% the contribution of the electroweak penguin 
operators thus making the central value of e'/e larger; 

• the uncertainty of the analysis is reduced as a consequence of the less effective 
cancelation between the gluon and electroweak penguin contributions; 

• the NLO O(p^) corrections are well under control and amount to a 10% effect on 
the predicted e'/e . 

Such an encouraging outcome makes us look forward to the next generation of ex¬ 
periments, under way at CERN, FNAL and DA<hNE, that will reach the sensitiveness 
of (1 -h 2) X 10“"^ and hopefully resolve the present uncertainty of the results of CERN 
(NA31) [8] 

Re {e'/e) = (23 ± 7) x 10"^ (1.1) 
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and FNAL (E731) [9] 


Re {e'/e) = (7.4 ± 6.0) x 10"^ . 


( 1 . 2 ) 


In the following subsections, before going into the detailed analysis, we summarize the 
notation and the framework used in our previous work. 


1.1 Basic Formulae 

The AS* = 1 quark effective lagrangian at a scale fi < rric can be written as [10] 


'^A5=l (h) T T|/j(p.) Qi{fi) 


(1.3) 


where i = I ,..., 10. The functions Zi(n) and are the Wilson coefficients (known to 
the NLO in the strong coupling expansion [11]) and Vij the Kobayashi-Maskawa (KM) 
matrix elements; r = -VtsV^jV^dV*^. 

The Qi are the effective four-quark operators obtained by integrating out in the stan¬ 
dard model the vector bosons and the heavy quarks t, b and c. They are by now the 
standard basis for the description of AS = 1 transitions and can be written as 


Qi = 

('Sa'W/9)v-A (^/3'^a)v—A ’ 

Q 2 = 

1 

> 

1 

> 

^3,5 = 

{sd)y_j^J2q (Q'5')v=fA ’ 

Q4,6 = 

{Sadp)y_^ X]g(Q'/35'a)V= fA j 

Q7,9 = 

2 ('S'^)v-A Sg (qq)v±a ’ 

Qspo = 

2 (^a(^0)v—A ^q 


(1.4) 


where a, f3 denote color indices (a, /5 = 1,..., W) and are quark charges. Color indices 
for the color singlet operators are omitted. The subscripts (V ± A) refer to 7^(1 ± 75 ) in 
the quark currents. 

Starting from the quark effective operators in eq. ( 1 . 4 ) and based on the xQM ap¬ 
proach we have given in ref. [ 12 ] a discussion of their bosonic representation and the 
determination of the corresponding 0{p^) chiral coefficients. In the present paper we 
complete the O(p^) part of the analysis by including a term proportional to the current 
quark masses previously neglected. 

Two additional operators enter when considering O(p^) matrix elements: 

Qn = (1 + 75)(1-75)] O'-Gd, ( 1 . 5 ) 

Qi 2 = [rridil + 'y^)+ms{l-'y5)] O'-F d. ( 1 . 6 ) 

The bosonization of the chromomagnetic operator Qn is discussed in ref. [ 6 ] where it 
is shown that the hrst non-vanishing contribution arises at 0{p‘^) and it is found to be 
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further suppressed by a factor m^/m\ compared to the naive expectation. Our complete 
NLO analysis shows the subleading role of these operators for tttt. 

For completeness we mention the so called self-penguin contribution to e'/e considered 
in ref. [13]. This contribution has not been included in any systematic NLO short-distance 
analysis. Since the effect on s'/e is of order 10“^, and thereby well within the errors 
quoted, it is omitted altogether. 

The quantity e'/e can be written as 


- = e^F2-.5o+^/4) Gpuj 


IIo-112 

UJ 


E 2 161 Re Aq 

where, referring to the AS = 1 quark lagrangian of eq. (1.3), 

1 


Ro = 

n2 = 


cos(5o 

1 

COs52 


'y ', yi RG(Qi)o (1 ^rj+ri' 

i 

Re((5i)2 , 


and 




(1.7) 


( 1 . 8 ) 

(1.9) 


( 1 . 10 ) 


Since Im = 0 according to the standard conventions, the short-distance component 
of s'/e is determined by the Wilson coefficients Ui. Following the approach of ref. [11], 
l/i(h) = 1/2 (h) = 0- As a consequence, the matrix elements of (^ 1,2 do not directly enter 
the determination of s'/e. 

Experimentally the phases 5i are obtained in terms of the tt-tt S-wave scattering 
lenght [14] at the rriK scale. The values so derived give to a few degrees uncertainty 
(5o — 37° and 82 — —8°, thus obtaining with good accuracy 


cos(5o ~ 0.8 

C0S(52 ~ 1.0 . (1-11) 

As a consequence the 1 = 0 amplitude includes a 20% enhancement from the rescattering 
phase. In addition, one has % — 52 — 7r/4 ~ 0 thus making the s'/e ratio approximately 
real. 

We take as input values 

~ 349 GeV■^ cu = 1/22.2 , = 0.25 ± 0.05 . (1.12) 

2|e|Re2lo 

The large value in eq. (1.12) for 1/uj comes from the AJ = 1/2 selection rule. In 
refs. [7, 15] we have shown that such a rule is well reproduced by the chiral quark model 
(xQM) evaluation of the hadronic matrix elements. 

The quantity includes the effect on 112 of the isospin-breaking mixing between 

7r° and the etas (see refs. [16] and [17]). Being this effect proportional to Ho, it is written 
for notational convenience as a part of the Ho contribution. Its uncertainty affects the 
error in the hnal estimate of s' je by about 10%. 
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As the precise values of Re Aq and oo depend on the choice of the xQM parameters— 
the quark and gluon condensates and the constituent quark mass M —we have used in 
the e'!e estimate the range of values that reproduce at 0{jA) the AJ = 1/2 rule with a 
20% accuracy [7]. 


1.2 The Chiral Quark Model 

Effective quark models of QCD can be derived in the framework of the extended Nambu- 
Jona-Lasinio (ENJL) model of chiral symmetry breaking [5]. Among them we have uti¬ 
lized the xQM [4], which is the mean-held approximation of the full ENJL and in which 
an effective interaction term 

is added to a QCD lagrangian whose propagating helds are the u, d, s quarks in the hxed 
backgronnd of soft gluons. 

In the factorization approximation, matrix elements of the four quark operators are 
written in terms of better known qnantities like qnark currents and densities. As a typical 
example, a matrix element of Qq can be written as 


\Qq\K^) = 2 (tt |u75(i|0)(7r^|sM|iL°) — 2 (7r’'"7r \dd\0){0\s''y5d\K^) 
(0|ss|0) — (0|(i(i|0)l {7r~^7r~\s'y5d\K^). 


(1.14) 


The matrix elements (building blocks) like (0| S 75 U \K~^{k)), {7r~^{p+)\sd\K~^{k)) and 
( 0 | 57 ^ (1 — 75 ) u |iL+(A;)) are then evalnated to O(p^) within the model and the four- 
qnark K tttt matrix elements are obtained by assembling the building blocks at the 
given order in momenta. 

Non-pertnrbative gluonic corrections, are included when calculating bnilding block 
matrix elements involving the SU{3) color generators T“: 


{0|st'‘T“(1-75)« |A'+(«:)) . 


(1.15) 


The use of quark propagators in gluon held background and the relation 

2 

= ^{-GG) (1.16) 

make it possible to express the non-factorizable gluonic contribntions in terms of the gln- 
onic vacnum condensate. The model thus parametrizes all current and density amplitudes 
in terms of {qq) , {asGG/ii) and M. The latter is interpreted as the constituent quark 
mass in mesonic helds. 

In order to restrict the values of the input parameters M, {qq) and {agGG/n) we 
have studied the AJ = 1/2 selection rule for non-leptonic kaon decay within the xQM, 
at 0(p"^) in the chiral expansion [7]. 
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The requirement of stability of the calculated isospin 0 and 2 amplitudes, togheter 
with minimal 75 -scheme dependence, selects the following range for M 

M = 200±20MeV, (1.17) 

which is in good agreement with the range found by htting radiative kaon decays [18]. 

For in the range 300 -7 380 MeV the O(p^) ht of ^ to their 

experimental values gives the ranges [7] 

{asGG/TT) = (334 ±4 MeV)^ , (1.18) 

(qq) = - (240 tlo MeV)^ (1.19) 

and 

M = 200 MeV , (1.20) 

for which the selection rule is reproduced with a 20 % accuracy. 

2 The AS = 1 Chiral Lagrangian 

One can write the 0{p^) AS = 1 chiral lagrangian in a form which makes the relation 
with the effective quark lagrangian of eq. (1.3) more transparent [12]: 

GW(g 7 , 8 )Tr (A^StAjs) 

+ G^"^\QyA [Tr (A^StA^SMts)+Tr (A^SA^StMSt)' 

+ G 8 (g 3 -io) Tr (A^D^St/l^s) 

+ GIAQi, 2 , 9 ,w) Tr (A?StZl^s) Tr (A^StD^s) 

+ G'^L(gi, 2 , 9 ,io) Tr (A^StD^s) Tr (A}StZl^s) 

+ GUiQ-rA Tr (A^D^SAjD^St) 

+ G’lj.iQ-rA Tr (A^St/l^s) Tr (AjED^St) 

+ Glj,{Q-rA [Tr (A?s) Tr (A'D^St/l^S S^) 

+ Tr (A^Zl^SZl^St S) Tr (A^^t)] , (2.1) 

where A® are combinations of Gell-Mann SU{3) matrices dehned by = AAjk 

and Af = diag(mti,mrf,ms). The matrix S is dehned as exp([jn(a;)^ and n(a;) = 
X]a-^“'^“(2^)/2, (a = 1,...,8) contains the pseudoscalar octet helds. The coupling / is 
identihed at LO with the octet decay constant. The covariant derivatives in eq. (2.1) are 
taken with respect to the external gauge helds whenever they are present. Other terms 
are possible, but they can be reduced to these by means of trace identities. 

The chiral coefficients G{Qi) can be determined by matching chiral amplitudes with 
those obtained in the xQM, via integration of the quark helds at the constituent quark 
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mass scale M. Given the 0{p^) chiral lagrangian we can fully compute the O(p^) meson- 
loop contributions to the hadronic matrix elements of interest. In refs. [7, 12, 15] one can 
hnd a detailed description of the approach we have adopted and the results obtained. 

The bosonization of the electroweak operators Qj^s starts with a momentum-indepen- 
dent term of order 0(p°) in the chiral expansion, this is the reason behind their 

importance in the determination of the ratio e'je . While the momentum-dependent 
0 {p^) corrections proportional to have been included in our previous numerical 

analysis, the momentum-independent correction proportional to was not determined 
in ref. [12] because, although of O(p^), it vanishes in the chiral limit ruq 0 and it was 
erroneously thought to give a negligible contribution. 

By matching the xQM calculation with that of the chiral lagrangian, for a convenient 
mesonic transition, one hnds 

G<’”>(Q,) = -af{qq)/N, 

G<”''(Q8) = -3f{qq} (2.2) 

In ref. [19] it was showed that should be proportional to the strong lagrangian 

coefficient Lg. This can be verihed via a strong chiral lagrangian calculation after inclusion 

of the wave-function and coupling-constant renormalization, as discussed in the appendix 
of ref. [7]. 

The determination in eq. (2.2) yields the following corrections to the Qg hadronic 
matrix elements for the 0 and 2 isospin projections of the tttt transitions: 

(Q 8 )o^^ = 2^3-^ (m^ -h Srud + 2 mu) 

(ms + 3md + 2m„) . (2.3) 

The one-loop chiral corrections for the 77° — tttt amplitudes induced by this term 
(including wave-function renormalization) are given by 

f2 /2 

= -^^^G(”^)(<57,8)(m, + 3mrf + 2m„) (0.759+ 0.444 ^ +0.163In/i2)(2.4) 

x-2 /H 

4,8 (j^) = -^^^G(^HQ 7 , 8 )(ms + 3md + 2mu)(3.48 + 0.240 i + 1.33lnju^) (2.5) 

where fi is in units of GeV and the labels (00,-1—) refer to the decay into two neutral or 
charged pions respectively. All the relevant chiral corrections for the other terms in the 
AS = 1 chiral lagrangian may be found in ref. [12]. 

The contribution of the G*^™^ term to the bosonization of the electroweak operators 
(^ 7,8 has the opposite sign with respect to the leading constant term. The effect of it, 
combined with the inclusion of the G^^ term (G^’^ are numerically negligible), and the 
0{p^) part of the xQM wavefunction renormalization, is to reduce by 20% the overall 
contribution of the electroweak operators. This determines a sizeable increase of the 
predicted central value of e'/e . 
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3 Hadronic Matrix Elements 


In order to perform the present analysis, we have computed the matrix elements of the 
relevant operator to O(p^), the next order in the chiral expansion. Such a computation 
was deemed necessary in order to further reduce the uncertainty in the estimate. 

The total matrix elements have the form 


(Qi(h))/ 







where Qi are the operators in eq. (1.4), 


{Q.), = ((7nr),|Q,|A-») . 


(3.1) 


(3.2) 


and are the xQM wavefunction renormalizations, whose expressions to O(p^) can 
be found in ref. [7]. The matrix elements so computed are expanded to O(p^), discarding 
higher order terms. The functions af (p) represent the scale dependent meson-loop cor¬ 
rections, including the mesonic wavefunction renormalization. They are dehned as the 
isospin projections of the and corrections computed in ref. [12], properly 

rescaled by factors of /,r// in order to replace / in the NLO evaluation. The scale 

dependence of the NLO part of the matrix elements is a consequence of the perturbative 
scale dependence of the current quark masses which enter at this order. 

A complete list of the hadronic matrix elements to order O(p^), namely {Qi)f^, can 
be found in ref. [12], with the proviso of replacing all occurences of /.„- by /. The NLO 
1 = 0 and 2 contributions to the tttt matrix elements are given by: 




-X 


3 

-ix 

3 




-1 + 


X. 


P - J^^{GG)(^G 




-2 + 


N, 




2 
A 

<^(GG> 

N, 

5 


/9g 


c 

(GG) 

X, 


{(3g + 37g) 


X, 


-X 


P - S{GG) {Pg - 7g) 


\NLO 


{Q 2 Y, 

—X 7 

AT / ’ 


- {QiX^ + (Qs) 


NLO 

0 


X, 
2X 
3 


, 5{GG) 

7 + ^:^7g 


X, 

.NLO 




ed) 


2 i^u Crf) 


1 

2 


V2XfJ + 2V6!jpjE 






(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 



(Qio)r° 

(Qll)r^ 

where e^ = 2/3, e^ = 


2^d{Q6)o^^ + 2 ^ + '^(CG) (^G + 37g) 

- 2^3 (e„ - Crf) {qq)'yE , 


' 2 ( 6 ^. erf) 


V^X r 


P + ^{GG)Pg 


Vg {cu - Crf) {qq)'yE , 


■^ed{Q3)Q 


NLO 


+ -{eu 


ed) (Qi) 


NLO 

0 


2 (^ 1)2 


NLO 


— 7,^d{Qi)N^ + o (e^ - erf) (< 22)0 


ATLO 


= - e„ 


ed) (Q 2 ) 


AfLO 

2 


T ^ -<0 

4 mj 




-1/3 are the up-, down-quark charges respectively, 
A' = V3/ (mjf - ml) 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


and the gluon-condensate correction (5 (gg) is given by 


TV, (asGG/n) 

2 167r2/^ 


The quantities (3, /5g, 7 , 7g and 7 ^ are dimensionless fnnctions of the mass parameters: 



(9 


/9g 

7 


7g 


m 


-h 2 m^ 


M 

IG -^ - 3^ K + 3m) 


m. 


m 


M 


1-6 



m; 


(m^ 


■ + 


m|- -|- 2m 


m; 

5ms + 17m (m^ — m)m 


m 

M 


1 - 12 


6 M 2 


(gg) 

P 


2MAI 

X, 


12M 12M {m\ — mp) 

m]^{2ms + 6m) — mpmg + 7m) 



2 m{ms 


m 


M{m]^ — ml 

+fl 


u 


{m\ — m 2 )A 
M 2 


m|- m^ 



4M2 


+ 


(m^ -f 5m) 2m{ms — m) 


2M 


mj^ — m^ 


3M/+ 

■ 

m| 


m|- m\{ms -l- 3m) — 2mpms -f m) 


2M 

ms —m 


(m 


K 


mi 


m\ — m2 


6M 


(3.22) 

(3.23) 


(3.24) 

(3.25) 
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lE 


{m\ + mlf + 3m^ 
4/MA2 




{m \ + ml) 

fK 


m- 


X2m\ + 13m^^ 


/A^ 


1 

7m 

fU 

{qq) 



(3.26) 


Concerning the gluonic operator Qu, we recall that, due to its chiral trasformation 
properties, it contributes only to the isospin zero amplitude and that {Qn)^‘^ = 0. The 
NLO matrix element of Qu is small compared to due to the additional suppres¬ 

sion factor m^/m\ [6]. 

We notice that strictly speaking at O(p^) also the meson two-loop corrections to the 
leading 0(p°) term of the operators (^ 7,8 should be included. However, these are chirally 
suppressed by 0(M^/A^) in comparison to the contributions computed in the xQM 

and will be neglected. 


3.1 The Bi Factors 


The effective factors 


^(0,2) ^ 


Re 




/0,2 




/ 0,2 


(3.27) 


give the ratios between our hadronic matrix elements and those of the VSA. They are a 
useful way of comparing different evaluations. In Table 1, we collect the factors for the 
relevant operators. Their values depend on the scale at which the matrix elements are 
evaluated an on the values of the xQM input parameters; moreover, they depend on the 
75 -scheme employed^. We have given in Table 1 a representative example of their values 
and variations in the t’ Hooft-Veltman (HV) scheme for 75 —that we employ throughout 
the analysis. 

In Table 1 we have indicated by “LO” the complete OijP') results, for / = f^, includ¬ 
ing the xQM gluonic corrections and the mass term proportional to (discussed in 
section 2 ), by “- 1 - y-loops” the sum of LO and chiral loop contributions, taking / at its 
renormalized value of 86 MeV (for a discussion on its determination and renormalization 
scheme dependence see ref. [7]), and by “-f NLO” the complete result which includes the 
effect of the O(p^) xQM corrections to the matrix elements here computed. 

The large values of Bf'l compared to B^l reflect the AJ = 1/2 rule in tttt 

decays which is well reproduced in our approach [7]. The evaluation of the penguin 
matrix elements {Q 3 ) and {Q 4 ) in the xQM leads to rather large Bi factors. In the case 
of Q 3 , the xQM result has the opposite sign of the VSA result and B^ is negative. This 
is the effect of the large non-factorizable gluonic corrections. 

The linear dependence in (qq) of the gluon penguin operators Qq and Qs in the xQM 
matrix elements (to be contrasted to the quadratic dependence of the VSA) is responsible 

^The 75 -scheine dependence enters at 0{p^) and has been discussed in ref. [3]. 
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LO 

-|- y-loops 

+ NLO 


4.2 

7.8 

9.5 

w 

1.3 

2.4 

2.9 


0.60 

0.58 

0.41 

Bs 

-0.62 

-1.9 

-2.3 

Ba 

1.0 

1.6 

1.9 

fig ~ Bq 

1.2 0.72 

2.0 H-1.2 

1.9 H-1.2 


0.91^-0.81 

2.6 H-2.4 

2.6 H-2.4 


1.8 

3.0 

3.6 


1.8 

3.7 

4.4 


0.81^-0.82 

0.95^-0.87 

0.94^-0.91 


0.60 

0.58 

0.41 


Table 1: The Bi factors in the xQM, leading order (LO), including meson-loop renormalizations 
(x-loops) and O(p^) terms (NLO). We have taken the gluon condensate at the central value of 
eq. (1.18), while the range given for B^-% corresponds to varying the quark condensate according 
to eq. (1.19). The results shown are given in the HV scheme for M = 200 MeV and / = 86 
MeV, except for the first column where f = fn has been used. 


for the sensitivity of q to the value of the quark condensate. 

The LO values for the electroweak penguin operators Qj^s are smaller than one because 
of the suppression induced by the O(p^) corrections to the leading 0(p°) term discussed 
in section 2. 

It is interesting to compare these results with those of lattice QCD. The lattice es¬ 
timate at the scale /r = 2 GeV for these operators gives = Bq = 1.0 ± 0.2 [20], 
= 0.58 ± 0.02, Bf'> = 0.81 ± 0.03 [21] and Bf^ = 0.62 ± 0.10 [20]. Even though 
the scales are different, these factors have been shown to depend very little on the energy 
scale [22]. We can therefore notice the qualitative agreement, at least for Bjg^g, on the 
fact that all values are smaller than the corresponding VSA. A lattice estimate of B^ 
would be a desirable check of the modeling of non-factorizable contributions in the xQM. 


4 The Mixing Parameter ImA^ 

In this section we update the determination of the combination of KM entries Im 
which is crucial in the study of e'/e . Such an update is required in the light of our recent 
O(p^) estimate of the parameter B^ given in ref. [7]. 

A range for Im Xt is determined from the experimental value of e as a function of 
rrit and the other relevant parameters involved in the theoretical estimate. We will use 
the recent NLO results for the QCD correction factors [23] and vary the AS = 2 
hadronic parameter Bx around the central value obtained within the xQM [7]. 
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In order to restrict the allowed values of Im A* we solve the equations 


£*^( 77 , p, \Vcb\, \ Vus\, AQCD,rnt,mc, Bk) = e 


(4.1) 


and 


2 I 2 
V +P = 


1 IK 


ub\ 


IK. 


+ 


l-p 1 


IK, 


IKK IKK 

|2 


1 IK 


td\ 




(4.2) 

(4.3) 


in terms of the rj and p parameters of the Wolfenstein parametrization of the KM matrix. 
We thus hnd the allowed region in the p — p plane, given rrit, rric and [24] 


e 

= (2.266 ± 0.023) X 10“^ 

(4,4) 

IK.I 

= 0.2205 ±0.0018 

(4.5) 

IKI 

= 0.040 ±0.003 

(4.6) 

IKKIKI 

= 0.08 ±0.02. 

(4.7) 


For |Kd| we use the bounds provided by the measured B^-B^ mixing according to the 
relation [1] 


\yu\ 


200 MeV 

'170 GeV' 

0.76 

[ AMb, 1 

0.5 

/O .55 


_ mt{mt) _ 


_0.50/ps_ 


/ Vb 


(4.8) 


with the ranges for the various parameters given in appendix A. 

The theoretical uncertainty on the hadronic matrix element controls a large 

part of the uncertainty on the determination of ImAj. For the renormalization group 
invariant parameter Bk we take the range found to O(p^) in the xQM [7]: 


Bk = 1.1 ± 0 . 2 , 


(4.9) 


where the error of the above direct determination includes also the whole range of values 
extracted independently from the study of the K^-Kg mass difference [7, 25]. 

The NLO order QCD corrections pi, 2,3 in the AS = 2 Wilson coefficient are computed 
by taking = 340 ± 40 MeV [24], mb(rnb) = 4.4, mc{mc) = 1.4 and = 175 ± 6 

GeV [26], which in LO corresponds to rritimw) = 177±7 GeV, where the running masses 
are given in the MS scheme. As an example, for central values of the input parameters 
we hnd at p = me 

Pi = 1.33 , p2 = 0.51 , p3 = 0.44. (4.10) 

Figs. 1 and 2 give the results of our analysis for the two extreme values of mt- the area 
enclosed by the two black circumferences represents the constraint of eq. (4.2), the area 
between the two gray (dashed) circumferences is allowed by the bounds from eq. (4.3); 
the area enclosed by the two solid parabolic curves represents the solution of eq. (4.1) 
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Figure 1: Range of allowed {r],p) values for 


169 GeV. See the text for explanations. 



Figure 2: Same as in Fig. 1 for = 181 GeV. 

with the range of Bk given in eq. (4.9) (notice that the upper parabolic curve corresponds 
to the minimal value of Vd, and vice versa for the lower curve). 

The gray region within the intersection of the curves is the range actually allowed 
after the correlation in Vd between eq. (4.1) and eq. (4.3) is taken into account. A 
further correlation is present when computing Im Xt from rj in eq. (4.11). 

This procedure determines the allowed range for 

\mXt-il\Vus\\Vd?. (4.11) 

For a flat variation of the input parameters, including Aqqj-,, we hnd 

0.62 X 10“^ < ImAt < 1.4 X 10“^ . (4.12) 

For comparison, values of Im Xt between 0.86 x lO”'^ and 1.7 x 10“^ are found in the most 
recent update of the Munich group [1] using Bk = 0.75 ± 0.15. The larger values of Bk 
we obtain reduce substantially the maximum values allowed for Im Xt. 
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5 Estimating e’je 


e7e 

2.5 
2 

1.5 
1 

0.5 
0 
-0.5 


Figure 3: Anatomy of e'je in units of 10“^ for central values of the input parameters: LO 
calculation (black), LO with chiral loops (half-tone), complete NLO result (gray). In the LO 
case we have taken f = fn, whereas in the remaining hystograms the renormalized value / = 86 
MeV has been used. 

We can now discuss our results for e'/e . Fig. 3 shows the impact on the hnal value 
of e'je of each operator for the representative central values of all input parameters 
and ImAt = 1.0 x 10“^. The same hgure also shows how the results vary from the 
LO predictions (black), once chiral loops (half-tone) and NLO (gray) corrections are 
included. The typical size of the xQM NLO corrections to the matrix elements of the 
quark operators is less than 10%. 

The contribution of the chromomagnetic penguin Qn is very small for two reasons, 
as discussed in ref. [6]: its matrix element is vanishing at the LO, and the NLO result 
turns out to be proportional to m^. Numerically, the matrix element of Qu is about 5% 
of the NLO corrections to the matrix element of Qq. 

The LO reduction of the contribution of the operator Qs caused by the O(p^) cor¬ 
rection terms is an important result of the present analysis. In fact, without them the 
contributions of Qq and Qs are at LO almost equal in absolute value thus leading to a 
vanishing value for s' je . With the inclusion of the 0{p^) corrections the hnal value of 
s' je turns out to be positive in the whole parameter range. It is important to remark 
that this result is also determined by the enhancement of the Qq matrix element due to 
the chiral loop corrections. In turn, this effect is related to the good ht of the amplitude 
AoiK^ —>■ tttt) which we obtain in the same framework [7]. 

In order to understand the uncertainty in our estimate, we now present a detailed 
study of the determination of e' je . The matrix elements depend on the values of the 
input parameters {qq) and (a^GG/Tr) , besides that of the constituent quark mass M 
characteristic of the xQM. 

The range of values for {qq) , (a^GG/vr) and M found in ref. [7] via the NLO best ht 
of the AJ = 1/2 rule are reported in eqs. (1.18)-(1.20). Since the value of (a^GG/vr) is 



Q3 Q4 Q5 06 Q7 Q8 Q9 QIO Qll Q12 All 


14 














Figure 4: e' je in units of 10 ^ as a function of {fiq) (GeV) and A = Im \t in units of 10 ^, 
for central values of the other input parameters. 



Figure 5: The dependence of e' je in units of 10“^ on the matching scale // as a function of 
(g(?) (GeV) and A = Im Aj in units of 10“^. The gray surface corresponds to /r = 0.8 GeV 
and the dark surface to ^ = 1 GeV. 

not very important in the physics of penguin operators we can safely consider the central 
value of eq. (1.18) as a hxed input for our numerical analysis of s' je . 

A relevant parameter in determining the size of s' je is ImA^, as discussed in the 
introduction. We use the overall range given by eq. (4.12). 

Fig. 4 shows s' je as a function of (gg) and Im A^ for central values of the other input 
parameters. The stability of the result can be gauged by comparing the variation of the 
surface in Fig. 4 as we vary 

• The matching scale /i, between 0.8 and 1 GeV: Fig. 5; 

• The values of Aqcd and mt in the ranges given in the appendix A: Fig. 6; 

• The value of M in the range of eq. (1.20): Fig. 7. 
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As we can see from Fig. 5, the dependence on the matching scale is very weak and 
vanishes within the range of {qq) shown. We consider the scale stability a snccess of onr 
approach. 



Figure 6: The dependence of e'/e (in units of 10 on Aqcd and mt as a function of {qq) 
(GeV) and A = Im \t in units of 10“^. The upper surface corresponds to Aq^j^ = 380 MeV and 
GeV while the lower surface corresponds to Aq^^ = 300 MeV and = 181 

GeV. 



Figure 7: The dependence of e'/e in units of 10“^ on the constituent quark mass M as a 
function of {qq) (GeV) and A = Im Xt in units of lO”'^. The upper surface corresponds to 
M = 197 MeV, Aq^j^ = 380 MeV and = 169 GeV, while the lower surface corresponds 

to M = 205 MeV, A^J.^ = 300 MeV and = 181 GeV. 

By varying Aqcd, (qq) and ImA^ in their respective ranges while keeping M hxed 
at its central value, we hnd (Fig. 6) 

e'/e = 1.7 X 10-3. (5.1) 

The hnal source of uncertainty arises from the value of M. If we include its variation 
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according to the range of eq. (1.20) we find (Fig. 7) 

e'/e = 1.7 X 10-^ (5.2) 

The result in eq. (5.2) represents the most conservative estimate of the error in our 
prediction. The dependence on the remaining input parameter in the NLO corrections 
is small and can be safely neglected. 

In comparing the present estimate with our previous one in ref. [3] (circa 1995), we 
may notice that the complete 0 {p^) analysis of the bosonization of the operator Qs has 
reduced the impact of the electroweak corrections and accordingly e'/e is never negative. 
In addition, the inclusion in the present analysis of the effects of the rescattering phases, 
eqs. (1.8)-(1.9), enhances the 1 = 0 channel, thus increasing the size of e'/e . As a 
consequence, the values of e'/e are now larger even though the NLO determination of the 
parameter Bk = 1.1 ± 0.2 has made the maximum value of Im At roughly 30% smaller. 

Finally, the updated new short-distance analysis, with the reduced range in Aqqq, 
makes the whole estimate more stable. 

5.1 Outlook 

Our analysis, based on the implementation of the xQM and chiral lagrangian methods, 
takes advantage of the observation that the AJ = 1/2 selection rule in kaon decays is 
well reproduced in terms of three basic parameters (the constituent quark mass M and 
the quark and gluon condensates) in terms of which all hadronic matrix elements of the 
AS* = 1 lagrangian can be expressed. We have used the best £t of the AJ = 1/2 selection 
rule in kaon decays to constrain the allowed ranges of M, (qq) and (GG) and have fed 
them in the analysis of e'/e based on the NLO O(p^) determination of all hadronic matrix 
elements. Values of e'/e positive and of the order of 10“^ are preferred, even though values 
of 0(10“^) cannot be excluded. 

In Fig. 8 we have summarized the present status of the theoretical predictions for 
e'/e. The two ranges we have reported for ref. [1] corresponds to taking two different 
determinations of m*, namely, from left to right, ms(1.3 GeV) = 150 ± 20 MeV and 
ms(1.3 GeV) = 100 ± 20 MeV. The varying of m* in the approach of ref. [1] is the 
analogue of varying {qq) in our LO matrix elements. 

The smaller uncertainty in the lattice result [2] is due to the authors’ Gaussian treat¬ 
ment of the uncertainties in the input parameters, to be contrasted to a flat scanning. 
The results of ref. [1] as well as ours can be made proportionally smaller by a similar 
treatment of the errors. 

Given the complexity of the computation, it is rewarding to find that so very different 
approaches lead to predictions that are reasonably in agreement with each other. 

We would like to conclude by briefly discussing three issues that arise in comparing 
our present analysis with that of refs. [1] and [2]: 

• The weight of the operator Q 4 , that is negligible in our estimate (Fig. 3), is made 
sizable in those of refs. [1] and [2] by the assumption that the coefficient B 4 can be 
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Figure 8: Present status of theoretical predictions and experimental values for s' je (in units 
of 10“^). The two horizontal bands bounded by long- and short-dashed lines show the \a 
experimental results of NA31 and E731 respectively. The vertical bars represent the theoretical 
predictions including the errors quoted by the authors. In gray we show our final estimate of 
eq. (5.2). 

as large as 6. Such a large value comes from considering the relation 

QA: = Q 2 — Q\ ^ Qz 1 (5-3) 

which is exact in the HV scheme, and the large value of the difference — (Qi)o 
that is required in order to account for the AJ = 1/2 rule. The xQM result 
of a negative Bj, shows that a large (< 22)0 ~ (Qi)o necessarily implies a 

correspondingly large value of i? 4 . A large value of makes s'/e smaller because 
the operator gives a contribution of the opposite sign with respect to that of Q%. 

• In ref. [1] only the terms proportional to and of the AS* = 1 chiral 

lagrangian are included in the matrix elements of the operators Qt^. While the 
effect of the missing terms is within the variation of the coefficient B^ that is 
there considered, the central value of s'/e may be affected by such a choice. 

• It is difficult to compare our approach to that of the lattice. We however notice that 
the contribution of the chiral term proportional to only starts at the level of 
the K ^ 271 amplitude and therefore is not included when computing the two-point 
K —>■ 71 transition, as currently done on the lattice. We would also like to stress 
the importance of always using the full O(p^) AS = 1 chiral lagrangian (2.1) in 
computing the K ^ 271 from the K ^ tt amplitude. 

• Finally, the 20% enhancement of the 1 = 0 amplitude due to the rescattering phase 
is included only in our analysis. 
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A Input Parameters 


parameter 

value 

sin^ 0 w{rnz) 

0.23 

mz 

91.187 GeV 

mw 

80.33 GeV 

(pole) 

ml ' 

175 ± 6 GeV 

mtinit) 

167 ± 6 GeV 

mt{mw) 

177 ± 7 GeV 

mbinib) 

4.4 GeV 

mcirric) 

1.4 GeV 

TTis (1 GeV) 

178 ± 18 MeV 

mu + md (1 GeV) 

12 ± 2.5 MeV 

A(4) 

^4qcd 

340 ± 40 MeV 

Vud 

0.9753 

Vus 

0.2205 ±0.0018 

\Vcb\ 

0.040 ±0.003 

\Vub/Vcb\ 

0.08 ±0.02 

VB 

0.55 ±0.01 


200 ± 40 MeV 

^Mb^ 

0464 ±0.018 ps-i 

Bk 

1.1 ±0.2 

Im At 

(1.0 ±0.4) X lO""* 

M 

200 MeV 

(m) 

-(240 I?)) MeV)3 

(asGG/TT) 

(334 ± 4 MeV)4 

f 

86 ± 13 MeV 

U = L+ 

92.4 MeV 

II 

+ 

113 MeV 

niTr = (m7r+ +m^o)/2 

138 MeV 

rriK = nifco 

498 MeV 

rur, 

548 MeV 

^•r)+rj’ 

0.25 ±0.05 

cos5q 

0.8 

cos 82 

1.0 


Table 2: Table of the numerical values of the input parameters. 
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